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Abstract 

We present an explicit numerical method to obtain the Cartan-Khaneja- 
Glaser decomposition of a general element G € SU(2 N ) in terms of its 
'Cartan' and 'non-Cartan' components. This effectively factors G in terms 
of group elements that belong in SU(2 n ) with n < N, a procedure that 
an be iterated down to n = 2. We show that every step reduces to 
solving the zeros of a matrix polynomial, obtained by truncation of the 
Baker-Campbell-Hausdorff formula, numerically. All computational tasks 
involved are straightforward and the overall truncation errors are well 
under control. 



1 Introduction 

In their seminal paper [K-G], Khaneja and Glaser presented a generic method 
to decompose 'large' unitary elements in terms of 'smaller' ones. The initial 
unitaries can be viewed as evolution operators of a multipartite system of spin- 
1/2's or as quantum algorithms acting on qubits. Experimentally it is only 
possible to control the interactions between a small number of subsystems at 
a time. Hence, this decomposition is of particular interest. It allows us to 
address questions such as how to optimize a dynamical evolution in terms of 
control procedures applied to a small number of spins or how a certain quantum 
algorithm can be obtained with the smallest possible number of experimentally 
feasible one and two qubit gates. 

In particular, Khaneja and Glaser show that any element of the Lie group 
G = SU(2 N ) is, up to local unitaries in SU(2)® N , determined by components 
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generated from certain Abclian subalgebras f)„ and f„, n = 2, ...,N, of the Lie 
algebra su(2 w ). This is achieved by employing iteratively the Cartan decompo- 
sition 

G = K'HK", 

where H is generated by ()„ (f„) and the factors K' and K" belong to the 
subgroup K C G generated by a particular subalgebra orthogonal to f)„ (f„). 
These relevant substructures are specified in terms of a fortunate choice of basis 
for su(2") that can be obtained recurrently for n = 2, N. 

The beauty of this result and its promising applications in quantum al- 
gorithms [V-W, S-B-M], control theory, quantum error correction [S] or the 
quantification of entanglement in multi-partite systems [O-S] have motivated 
the search for a constructive method to perform the decompositions for any 
given matrix G E SU(2 N ). Although substantial work has been done on the 
first non-trivial instance SU(4) [K-C, H-V-C], little seems to be known so far 
for the higher N case [Bu]. 

Here we address the decomposition problem for the general group SU(2 N ). 
Employing a convenient truncation of the Baker-Campbcll-Hausdorff (BCH) 
formula, we show that the problem allows for a numerical algorithm to calculate 
all such KHK decompositions with controlled error. Hence, we can perform 
the full Khaneja-Glaser decomposition of a general element G G SU(2 N ) with 
arbitrary computational precision. 

This article is organized as follows. In Section 2 we briefly review the 
Khaneja-Glaser decomposition [K-G] and establish the formalism for our ap- 
proach. In Section 3 we reinterpret the problem in terms of the BCH formula; 
we explain how a truncation of the BCH series renders our problem solvable by 
straightforward numerical tasks. In Section 4 we give a user-friendly summary 
of the essential steps involved and we finally conclude in Section 5. 



2 The Khaneja-Glaser decomposition 

We consider the compact semi-simple Lie group G = SU(2 N ) and a particular 
compact closed subgroup K C G; we denote by g = su(2 N ) the Lie algebra of G 
and by t C g the Lie algebra of K understood as a subalgebra of g. Since G is 
semi-simple, the Killing form (., .) is non-degenerate and defines a bi-invariant 
metric on G. Hence, we can define m = t to be the orthogonal complement of 
t with respect to the metric. Notice that, in general, the vector space m C g is 
not a subalgebra. Since m is determined by t and the Killing metric, we shall 
refer to this structure as the Lie algebra pair (fl, 6). 

We shall adopt the following typographic conventions in most cases 

G capital bold group or subgroup 

G capital element of G 

g German fraktur Lie algebra or subspace 

g normal element of g 
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The only exception to the choice of fonts just stated will be the familiar 
Pauli matrices, seen as elements of su(2), which will be denoted by majuscules 



2\10y 2 ^ i / 2 V -1 

and the 2x2 identity matrix I = ^ J 1 ) ' ^^ rLa ^' wc snau use an abbreviated 
normalized notation for their tensor products given by 

(2\ fe ~ 1 A = X Y Z 

t) A l ®A 2 ®...®B j ®...®A k , B l =/ ' ' 

For example, YXI = 2 ( ? ~q ) ( 1 ) ® ( 1 
2.1 Cartan decomposition 

Definition 1 ^4 Cartan decomposition of g is an orthogonal split of g 

g = t®m 

given by a Lie algebra pair (0,6) satisfying the commutation relations 

[t,f]ct, [m,«]cm, [m,m]cfi. (1) 
In this case (g, 6) is called a symmetric Lie algebra pair. 

Remark 2 The apparently artificial conditions in the above definition have an 
interpretation in Riemannian geometry: g = t © m is a Cartan decomposition 
if and only if the quotient manifold G /K = cxp (m) is a globally Riemannian 
symmetric space [K-N]. Such a space possesses a canonical global involution 
(i.e., an automorphism s of the space onto itself such that s 2 = I) which induces 
naturally a linear involution s* on g that preserves the Lie algebra structure [.,.]. 
Since s* squares to identity, its eigenvalues must be ±1, and the subspaces 6 and 
m emerge respectively as the +1 and —1 - eigenspaces o/s*. 

We can start exploring the Cartan decomposition g — 6 © m by noticing that 
if \) C m is a subalgebra of g then, by (1), f) is automatically Abelian. Since m 
itself is not in general a subalgebra of g, it is natural to look for a maximal 
(Abelian) subalgebra \) C m. 

Definition 3 A Cartan subalgebra of the pair (g, 6) is a maximal (inm) Abelian 
subalgebra f) C m. 

From now on we shall assume that the Cartan subalgebra f) refers to the 
pair (g, t) unless stated otherwise. The following Proposition shows that the 
whole m is obtained from t) by the adjoint action of K and that there is only 
one such f) up to this action. In the context of the particular application we 
have in mind, this means \) carries the 'essential' information about m. 
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Proposition 4 Let f) and f)' be two Cartan subalgebras; then 
1. m= |J Ad/f ([)); 

JT6K 

£ t)' = Ad K (I)) for some element K S K. 

Proof. Please see Appendix B. ■ 

Denote H = exp(f)) C G the subgroup generated by f). The Cartan decom- 
position theorem states that any group element G £ G can be written as an 
element H e H together with left- and right-multiplications by elements of K: 

Theorem 5 (Cartan decomposition) The Cartan decomposition g = t © m 
induces a Cartan decomposition of the group G, 

G = KHK, (2) 

where H = cxp([)). 

The Cartan decomposition of a given element G € G has the form 

G = KqK^HKI = K M, (3) 

M 

where K ,Ki G K, H <G H and M = KiHK\ e exp (m). 

Proof. Since G/K = exp (m), there exist K e K cmd M e exp (m) such that 
G = K M. Let m — log (M) e m; from Proposition 4, item I, there exists 
Ki <E K such that Ad/^ m = h e f), so 

G = i^o ex P (™) = -^0 exp ^Ad^t hj = K KiCxp (h)K\. 

H 



2.2 The Khaneja-Glaser basis 

The Khaneja-Glaser basis [K-G] for arbitrary su (2™) makes explicit all the struc- 
tures which concern us. In particular, the splitting g = 6 © m and the Cartan 
subalgebras f)„ and f„ are manifest. 
We start with a familiar example: 

Example 6 For n = 2, we propose basis elements for the Lie algebra su (4) and 
organize them into subspaces m.2 and £2 as follows: 
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su (4) 


m 2 

XX XY XZ 
YY YZ YX 
ZZ ZX ZY 

u 

f)2 


h 

XI IX 
YI IY 
ZI IZ 



f) 2 = span {XX,YY,ZZ} 
h = {0} 

Now we obtain the Khaneja-Glaser basis for su(2") by a relatively simple 
iteration process starting from su (4), as summarized in the following diagram: 



su(2 n ) 




I n ~ l ®X I n ~ l ®Y 
su(2"- 1 )®X su(2"- 1 )®y 
U 


jn- 






su (2 n - 1 )®Z 
U 

fn 


fn 


(4 

su(2™- 1 ) ®I 

fin.O 



jn-1 = j®(n-l) = 

n-1 

NB.: Terms of the form su (2™ _1 ) 0^4 denote the set obtained by tensoring each 
element of su(2™ _1 ) with the matrix A = X,Y,X,I; in all cases we consider 
the span of the basis elements over M. 

Notice that for n > 2 two successive Cartan decompositions can be per- 
formed. First, the expected one referring to the pair (su(2™) ,{*„), with Cartan 
subalgebra f)„. In the terms of Theorem 5, this means we can write G <E SU (2 n ) 
as G = K'HK" with H e cxp (f)„) and K',K" G exp(t„). However, the de- 
composition of 6„ given by diagram (4) is 

6„ = fin.i 9 6„,o 9 span {r^Z} ~ su (2"- 1 ) © su (2"- 1 ) © u (1) , 

where both t n ,i and 6 ni o are canonically isomorphic to su(2™ _1 ). Since our 
aim is to iteratively decompose the factors K^) until they eventually reduce to 
'local' unitarics in SU (2)® N and non-local 'Cartan' factors, we would expect 
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t n to generate something of the form SU (2™ x ) ® SU (2). Thus, there is a 
su(2 Tl_1 ) component too many in 6„ that we need to factor away in order to 
define the complete recurrence step. 

A closer look at fi„ reveals another Lie Algebra pair, thereby clearing the way 
for a second Cartan decomposition: we just have to leave aside the 'complex 
phase' generated by I n ~ 1 Z (see diagram (4)) that can be seen as a 'local' trans- 
formation under the inclusion U(l) ^ SU (2). Let t n = 6n,i © 6n,o denote the 
subalgebra obtained from t n in this manner, so that i n = £n©u (1) t n ®su (2). 
Accordingly, given a group clement K = cxp (fc) G K, let us write K — exp ^A;^ , 

where k € t n is obtained from k € i by eliminating the component spanned by 
I n ~ 1 Z. This is well defined as I n ~ 1 Z commutes with every element in t n . 

It is now straightforward to check that (t n ,t n ,o*J 1S al so a Lie algebra pair 
[K-G], whose Cartan subalgebra we call f„. Hence, we can apply Theorem 5 
again to decompose the factors K' and K" into elements of exp (f„) together 
with left and right multiplication by some new factors generated by £ n ,o- This 
time the orthogonal subalgebra t n .o = su(2™ _1 ) ® I is canonically isomorphic 
to su(2"- 1 ), so it generates exp (t nfi @ span {r^Z}) ~ SU{2 n ' 1 ) ®SU{2). 
Thus, we have accomplished the the complete nth recurrence step that yields 
the decomposition 

g = kMfMkV>hkV>fV>kV\ 

with Ftf) e cxp (f„) , H e cxp (f)„) and e SU (2™- 1 ) ® SU (2). 

Note, finally, that we are particularly interested in the 'Cartan' factors, 
i.e. those generated by the Cartan subalgebras f)„ of (su (2™) , t n ) and f n of 

(^n^nfl^ji that emerge in each step. It is thus convenient to know explicitly a 
set of basis elements for each of these subalgebras. This can be achieved by the 
following recurrence formula, starting from = span {XX, YY, ZZ}, 



f) n = span a (n) 
f n = span b (n) 



n = 2,...,N 



(5) 



a(2) = {XX,YY,ZZ},b (2) = {0} 

n 

s (n) = (J a (j) ® f-J 

a(n + l) = {/ n ,s(n)}8)X 
b(n + l) = {s(n)} ®Z 



The Example below illustrates all the above constructions for the first non- 
trivial case su (8): 
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Example 7 n = 3 : su (8) 











su (8) 




























h 






IIX 






IIY 








IIZ 




su (4) <g> X 


su (4) <g> y 


su 


(4)8 


z 


t 


', 


su 


(4)6 
















«3,1 










^3,0 




U 

[)3 












U 
fs 













f) 3 = span {IIX, XXX, YYX,ZZX} 
f 3 = span{XXZ,rrZ,ZZZ} 

2.3 The Baker-Campbell-Hausdorff formula 

The matrix Lie algebra q is noncommutative and thus, for general elements 
a,b E g, the product of exponentials e a e b does not coincide with the exponential 
of their sum, e a+b . In fact the expression for log (e a e b ) has an infinite series 
of correction terms and is known as the BCH formula, after Baker-Campbell- 
Hausdorff. 

Although the original formula was rather complicated and computationally 
unpractical, a remarkable simplification made by Dynkin [D, B] expresses all 
the terms as successive commutators of a and b 

oo 

io g ( e a e ")= X^M). 

Here Tjj (a, b) denotes the homogeneous term of degree i in a and degree j in 
b; its expression is 

T lJ (a,b) = -^— Y \a il b j \..a ik b jk ] , 

' 3K ' i + J .. . , k n\jil..i k \j k \ L J 

•■•,»* ,3k) 

where we abbreviate [a n 6 J1 ...a lfc 6 :,fc ] = [a, ...[a,[b, ...[b,[a, ...[b,[a, ...[a,[b, ...,b]...] 

il jl ik jk 

and the sum ranges over all possible 2fc-uples of non-negative integers (ii, ji, ik, 
such that 

k k 



~Y^ic = i, ^2jc=j and i c + j c > 0. 



c=l c=l 
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The first few terms are 

log (e a e h ) = a + b + \ [a, b] + ± [a, [a, b}} + ± [b, [b, a}} + 

+ ±[a,[b,[a,b]]] + m(~) (6) 

and the higher order coefheients after decrease quickly (see e.g. [R]). This 
will motivate us later on to perform convenient truncations on this convergent 
scries. 



3 Numerical algorithm for the KHK decompo- 
sition 

In this Section we will develop a technique that allows the explicit numerical 
calculation of the components of a general group element G G G under Cartan 
decomposition. The idea is to consider the Cartan decomposition (3) in the 
light of the BCH expansion (6). Let g G q, m G m and k G fi be the generators 
of G, M and K , respectively. Then (3) reads 

G = e9 = e k e m . (7) 

Since the matrix G is given, (3) shows that k can be obtained from m (and 
vice-versa) by 

k = k{m) =log(Ge- m ). (8) 
Hence, the decomposition problem is reduced to finding m. 

3.1 Determining m 

First, taking logarithms in (3), we obtain 

.g = log(e fe e m ). (9) 

We then apply (6) to expand g in terms of successive brackets of k and m. 
In the light of Remark 2 we can easily deduce that each of the brackets belongs 
in either t or m. Hence, the expansion is split into two orthogonal components, 

g = k + m+ ^ [k,m] + [k, [k,m\] + [m, [m,k]] + [k, [m, [k,m]]} + ... 

S v ' S v ' S v ' S v ' 

em em ee ee 

= 9t + 9m, 



where 



9i = k+^[m,[m,k]} + ^[k,[m,[k,m}]} + ... e «, (10) 
9m = m+ - [k,m] + — [k, [k,m]} + ... e m. (11) 
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Note that computing g( and g m from g is a straightforward task since the 
Khaneja-Glaser basis (4) makes explicit the partition g = 6 © m. 

At this stage we can use (3) to eliminate k — k (to) in either of the equations 
(10,11). Choosing (11) whose first few terms are simpler we obtain 

9m = g m (m) = m + i [k(m),m] + [k (m) , [k(m),m\] + (12) 

which is an infinite series with rapidly decreasing coefficients. As g m (to) is a 
converging series we can truncate it so that the resulting equation will provide an 
approximation of to with an error that decreases by including higher commutator 
terms. If we call P p (m) the truncation that includes all terms with up to p 
commutators, i.e. 

ffm (m) — P P (m) + further terms, (13) 
we can in principle solve the equation 

P p (to) ee P p (to) -g m = (14) 

with respect to the single matrix variable to. However, the expression of k (to) 
given by (3) is rather complicated. So we propose using again the BCH expan- 
sion to obtain 

k (to) = log(e 9 e" m ) = g - m - ^ [g, to] - [g, [g, to]] + [to, [to, g]] H(15) 

As before, we can truncate (15) to a term that includes all q-th order brackets. 
This yields a polynomial Q q (to) that approximates k (to) as well as we desire 
at the cost of taking extra high-order commutators 

k (to) = Q q (to) + further terms. (16) 

After both truncations, equation (14) is approximated by 

= P p (to) ~ -g m + to + i [Q q (to) , to] + ^ [Q q (to) , [Q q (m) , to]] + 

where P p (m) is now a polynomial in one matrix variable with matrix coefficients. 
Our problem of finding to has thus been reduced to finding the zeros of a poly- 
nomial. This can be easily performed with a numerical algorithm [D-P-M]. The 
accuracy of the result can be increased by including more terms in the truncated 
series P p and Q q , i.e. increasing p and q. Specifically, we prove in Appendix A 
that, in this way, the errors in determining to by the truncation procedure are 
well under control. 

Example 8 Take the "first order" truncations p = 1, q = 1: 
Pi (m) = -g m + to + i [k (to) , to] , 
Qi(m) = g-m-^{g,m}. 
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Then, approximating k (to) ~ Qi (m), we obtain 



P 1 (to) 



m 



-.9m + m + ^ [g, to] - ^ [m, [to, 3]] . 



3.2 The M = K X EK[ decomposition 

Once to is known, it remains to find the subgroup element K\ e K whose 
adjoint action rotates M onto H = e h , h G f). 

Lemma 9 TTie following properties are associated to H = exp (f)): 

J. H is a iorws (compact connected Abelian Lie subgroup) of G; 

2. any vector v G f) whose l-parameter subgroup {exp (to)} is dense in H is 
centralized in m jwsi fry f) : 



{u G m I = 0} = t) 



(17) 



Proof, see e.g. [W, §8.6]. ■ 

As shown in Appendix B, the first necessary ingredient to perform the decom- 
position M — KiHK\ is some vector v € t) that generates a dense l-parameter 
subgroup exp (to) C f). This may seem abstract, but since we have an explicit 
basis (5) for f), it suffices to take any irrational combination of the Cartan gen- 
erators. 

Example 10 in.su (8), take e.g. 

v = IIX + TV. XXX + tt 2 .YYX + tv 3 .ZZX. 

The reluctant reader may verify that, indeed, the centralizer of such v in m is 
just f). 

Now, we may define f VjTn as in the Appendix B to be given by 

f v;m (K) = (v, Ad K (to)) = ^ C c ad Ctv a (K^mK) b (18) 

a.b,c,d 

and recover K\ numerically as a minimum of /„ lTO . 

We conclude that to can be rotated into h = Ad^i (to). Thus we have 
completed the decomposition 



.-. G = K Kie K{ 
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4 Step-by-step summary 

What we have described so far consists of the main building blocks necessary 
to perform the Khaneja-Glaser decomposition. Here, we will summarize all the 
steps one needs to take when given an arbitrary unitary G £ SU(2 n ). 

1. Calculate its (matrix) logarithm g = log (G) £ g — su (2 n ). 

2. Compute the Khaneja-Glaser basis following the recurrence in diagram 
(4); take g m , the component of g on the subspace m„. 

3. Truncate (12) including p-th commutators to get P p (to); let P p (to) = 
Pp(m) -g m . 

4. Truncate (15) including q-th commutator to get Q q (to), as in (16). 

5. Replace Q q (m) for k (to) in the expression of P p (to) obtained in 3, so 
that P p (to) becomes a polynomial in to. 

6. Solve the zeros of P p (to) to get a solution to to (14). 

7. Use to from item 6 above to calculate K = Ge~ m as in (3). 

8. Compute \) n following (5); order its elements {uj} e.g. alphabetically and 
define v = Y] n 3 Uj to satisfy the density hypothesis of Lemma 9. 

9. Use to and v to define f v , m (if) = ^ C c ad C^ c v a {K^mK) h as in (18); 

a,b,c,d 

minimize / on K = cxp (fi„) to find K\. 

10. Calculate h = K\mKi, and thus H — exp(h). 

11. Assembling the results from items 7, 9. and 10, obtain 

G = K Q KxHK\. 

12. Repeat the above steps for G = K Ki and then for G = Kl, replacing 

in — > 6n,0) ™n - * in,l an d fyn — > fn- 

13. Items 11 and 12 yield the decomposition 

g = kWfv>kv>hkv>fv>km, 

with e cxp (f„) , ff £ cxp ([)„) and JfCi) g SU (2™- 1 ) ® SI/ (2). 

14. Decrease n — > n — 1 and iterate this process to further decompose each 
factor if W £ St/ (2™- 1 ) ® SU (2) until they all reduce to a product of 

Cartan factors and H$ and local unitarics in SU (2) <g ' N . 

NB.: As far as accuracy in step 6 is concerned, tasks 3 to 5 should be 
performed in the light of Appendix A. Namely, truncations at higher order 
should be tried until numerical errors are satisfactory, which will happen after 
a finite number of attempts. 
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5 Conclusions 



As the advances in quantum technologies move beyond the control of one or two 
spins, or qubits, it is important to minimize the overall cost of processing quan- 
tum information. The Khaneja-Glaser decomposition of SU (2 N ) offers an upper 
bound for this optimization procedure given by 4 JV ~ 1 multi-local SU(2)® N ro- 
tations together with ^l^ -1 — 1 purely entangling operations. The latter can 
be reduced, if desired, into bipartite interactions, or two-qubit gates. More- 
over, in [N], Nielsen gives lower bounds for such optimization. Here we exploit 
the Khaneja-Glaser approach to build a constructive method for decomposing 
a general unitary in terms of its local unitary components. Abstract as it may 
seem, the decomposition problem can be cast in such a way that can be easily 
solved by a numerical algorithm that can be found at 
http://cam.qubit.org/users/iiannis/liesolvellJ.tar.gz. 

Finally, one should notice that the solutions we obtained are not necessarily 
unique. In general, neither the zeros of the matrix polynomials, P p (m), nor the 
minima of the functions, /, are unique. In particular, viewing the minimization 
procedure from the equivalent point of view of the diagonalization of the ma- 
trix to, where K\ is constructed out of the eigenvectors of m, there arc many 
equivalent solutions depending on the particular ordering of the eigenvectors. 
Moreover, one should also take into account that the exponential function has 
a natural 2tt periodicity and the adjoint action is Z 2 -symmetric [K-C, H-V-C]. 
While our approach is not concerned with the actual parametrizations of the 
group elements, this is an important issue which should be addressed in the 
future. 

Acknowledgements 

We would like to thank Gregorio Malajovich for very fruitful correspondence 
and for writing the program to find the zeros of our polynomials. We also thank 
Tim Pcrutz for drawing our attention to Rouche's theorem. This work was 
partially supported by the Royal Society. 



12 



A Appendix: Accuracy of BCH truncations 

Here, a generalised version of Rouche's theorem [L] is employed to show that 
the truncations performed in Subsection 3.1 yield a rigorous approximation for 
the zeros of (12). 

Theorem 11 Let ip,ip : C r — > C r be holomorphic functions and DcC be an 
open domain such that neither <p nor ip have zeros on dD; if 

\<p (m) — ip (m)\ < \<p (m)\ + \ip (m)\ , Vto e dD, (19) 

then (p and ip have the same number of zeros in D. 

NB.: We adopt, e.g., the norm \ip\ = maxjj \iPij\i but the argument holds 
for any L p -norm. 

Let r = 2 2n be the number of entries of a matrix to G su(2"), seen as 
a complex vector. Consider then a BCH-type series ip (to) and its truncated 
version tp (to) = P p (to) containing all its terms of up to p successive brackets, 
calling R p (m) the truncation remainder 

ip ( m ) = P p (to) + R p (to) . 

Suppose to G su (2™) is a zero of P p (m); then ip (to) will also have a zero inside 
the polydisc D — As(m) C C r of radius 8 > about rh if the following (stronger) 
instance of (19) holds 

\Rp{m)\<\P p {m)\, \fmedA 5 (m). (20) 

In other words, rh approximates at least one zero of ip (to) with error inferior to 
an arbitrarily chosen 8. 

All we have to show is that condition (20) holds for suitably large p; this is a 
relatively straightforward consequence of the uniform convergence of the BCH 
series tp (to) = lim P p (to), as we will now see. For any (small) £o > 0, there is 

p— »oo 

a po such that 

p > p a \P P (to) - P po (to)| < e , Vm, (21a) 

hence 

\P p (m)\>\P po (m)\-e Q , Vm. (21b) 

In particular, if to is a zero of P po (to), equation (21a) restricted to the bound- 
ary <9A,5(too) implies that all polynomials P p (to), p > p , also have at least one 
zero inside As (mo), by Rouche's theorem. 

On the other hand, convergence also implies lim \R p (m)\ = 0, hence, for a 

p— »oo 

given e\ > 0, there is p\ such that 

p > px => \Rp (m)\ < ei, Vm. (22) 
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Set £1 = min \P pa (in) \ — e ; assuming this is positive (if not, take a larger 

m£9A { (mo) 

Po for a smaller £o), take some p > max{po,Pi} and find mi G A^rno) such 
that P p (mi) — 0. Then consider any smaller polydisc Ag'(mi) C A^too), and 
restrict equations (216) and (22) to its boundary. We obtain 

\R p (m)\ <\P p (m)\, Vm e dAy(mi), 

thus, by Rouche's theorem again, <p (m) — P p (m) + R p (m) has a zero inside 
Ay(mi) c Aj(m ). 



B Appendix: Proof of Proposition 4 



In this Appendix we will develop the proof of Proposition 4, following [W, §8.3]. 
This argument contains most of the crucial elements to understanding the KHK 
decomposition in detail. 
Proof of Proposition 4- 

1 . Given m G m, we want to find K\ G K whose action rotates m onto some 
element h of the Cartan subalgcbra t). First, take any v G f) that generates 
a dense 1 -parameter subgroup exp (tv) G H = exp(f)), as in Lemma 9 
and define the function 



fv,: 



K 



/(*) 



(23) 



(w, Ad/f (m)) 



where (a, b) = tr (ad a adf,) is the Killing form on q. Since / is continuous 
and K is compact, / admits an absolute minimum Ki G K. If we consider 
a local perturbation of K~i by e tk , for any k G 6, we have 



= 



t=o 



(v, Ad e tfc Ad^i (to)) 



t=o 



k, Ad Kl (m) 



tr (ad„ ad[ fei/l ]) 



= tr (ad„ (ad fe ad/, - ad h adfe)) (c = t ' ) tr (ad [M] ad fe ) 
([ft,u],A;) = 0, Vfc G £. 

But the Killing form is non-degenerate on 6, so we must have [h, v] = 
and thus ft. G f), as w is centralized by f) (2); hence, we have shown that 

to = Ad K t (ft) , ft G t). 

2. By 1. there exists if G K such that 

Ad K (v) G t)'. 

Now take centralizcrs on both sides. 
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